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Abstract. We present a theoretical study of bose condensation of non-interacting bosons in finite lattices 
in quartic potentials in one, two, and three dimensions. We investigate dimensionality effects and quartic 
potential effects on single boson density of energy states, condensation temperature, condensate fraction, 
and specific heat. The results obtained are compared with corresponding results for lattice bosons in 
harmonic traps. 
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1 Introduction of the periodic lattice potential and an overall confining 

harmonic potential. Many theoretical studies [ WTUlfTTlfTSl 



O 

Over last few years, bosons and fermions in optical lattices 
O ■ HHHITllTnEZllIHllIllin] of such lattice bosons in har- 

. have emerged as important controllable systems for inves- 
• *h . monic confinement have appeared in recent years. Unlike 

rS \ tigations into several properties of quantum many-particle 

^ , the case of lattice bosons in harmonic confinement, the 

systems pQ. These are clean systems in which experimen- 

effect of anharmonic potentials on the properties of lat- 
tice bosons have begun to be explored only recently [2"T]. 
In that work Gygi and collaborators studied zero temper- 
ature properties of strongly interacting lattice bosons in 
a quartic trap. It has been suggested|21j that it is experi- 
mentally possible to create an optical lattice in a quartic 
trap employing a combination of red-detuned and blue- 
detuned Gaussian laser beams. Finite temperature prop- 



talists have achieved great control over a wide range 
of particle numbers, particle hopping, and strength and 
sign of inter-particle interactions. Experimental groups 
have conducted extensive studies [2,3,4,5,6,7,8^ of sev- 
eral properties of many-boson systems in one, two, and 
three dimensional optical lattices. The many-boson sys- 
tem in these experiments is under the combined influence 
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erties of lattice bosons in a quartic trap is of considerable tonian of the many-boson system we consider is 
interest in this context. 

H = -t^2 (4 c i + 4 Ci ) + 12 v - ^12 ni ' w 

In this paper, we present a theoretical study of lattice <ij> v ' i , t 

bosons in a quartic potential in one, two, and three di- ..... . 

where t is the kinetic energy gam when a boson hop from 

mensions (Id, 2d, and 3d). We consider non-interacting . . .... .... + 

site i to its nearest neighbor site j in the optical lattice, c\ 

bosons in a periodic lattice in Id in a Id quartic poten- 

is the boson creation operator, V(i) is the quartic poten- 
tial, a square lattice in a 2d quartic potential, and a cubic ... + 

tial at site t, rij = c]c, the boson number operator, and /i 

lattice in a 3d quartic potential. We study the effects of . 

the chemical potential. The forms of the quartic potentials 

used are: V(i) = Q xf in Id, V(i) — Q (xf + yf) in 2d, and 

V(i) = Q (xf + yf + zf) in 3d. We first obtain the matrix 



the potential on one-boson density of energy states (DOS) 

used are: V(i) Q .- 1 ■ 1 ! '.>—'>'■■ 
and the temperature dependence of ground state occu- TA/ . N ^ , 4 4 4 



pancy and specific- heat. We compare the results obtained . ... 

representation of the system Hamiltonian m a site basis. 

for lattice bosons in quartic traps to the corresponding . . 

We numerically diagonahze it to obtain energy levels of 

results for lattice bosons in harmonic trapsEDl. This work 

a lattice boson. We have used open boundary conditions. 

is presented in the Sections 2 and 3, and conclusions are . . . 

We have chosen lattice sizes large enough so that finite 

given in Section 4. In the work presented in the following , , 

size effects are absent m the results presented m the next 

sections, we have not included the effect of boson-boson in- . . . 

section. The lattice sizes were fixed by finding the lattice 

teraction (U). Our results would approximately also hold 

size beyond which results remain unchanged with further 

in the weak interaction regime (U << t, where t is the . . 

increase of lattice size. I his depends on the magnitude of 

boson hopping energy) where the interaction induced de- . ....... 

the quartic potential strength since it decides the spread 

pletion effects are not significant 22 . A weakly interacting 

of the boson distribution in the lattice for a given value 

regime may be achieved by adjusting the lattice potential . . 

of t and temperature. The energy levels (Ei) obtained for 

depth to a low value as has already been done for lattice , , , , . . . . . 

a boson m these large lattices are used m calculations of 

bosons in quadratic traps fUS]. , n , , . n , 

the DOS, ground state occupancy, and the specific heat. 

The chemical potential and boson populations in the var- 
ious energy levels are calculated using the boson num- 

2 Model and method 

ber equation: N = }_, i=0 A where Eq and k, m are 
the lowest and the highest single boson energy levels and 
In this Section, we give a brief presentation of the model of N(Ei) — \/[exp[(3 (Ei — /i)] — 1] in which (3 = l/fc^T with 
the system and the method followed in the calculations of ks the Boltzmann constant and T the temperature. The 
various properties presented in later sections. The Hamil- specific heat is calculated from the temperature derivative 
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of total energy (E tat — J2iLoN(Ei)Ei). All energies are a 2d square lattice in a 2d quartic potential (Fig. 2(a)) 
measured in units of t. and the corresponding results for the harmonic case(Fig. 

2(b)). In the quartic potential case, we find that the van- 
_ „ , ■ r^- ■ Hove singularity is strongly suppressed by the confinement 

3 Results and Discussion 6 J 6 J 

potential. In comparison, the Van Hove singularity is de- 
3.1 One-boson density of energy states stroyed by the harmonic potential [23]. In both cases the 



In this section, in Figs. 1-4, we present our results on 
the one-boson density of energy states (DOS) for a lat- 
tice bosons in a quartic trap and compare with the corre- 
sponding DOS for a lattice boson in a harmonic trap for 
which: V(i) = K(x 2 + yf + zf). In Fig. 1, we have exhib- 
ited the DOS for bosons in optical lattices with harmonic 
and quartic confinement potentials for several values of 
q = Qa A and k — Ka 2 , where a is lattice constant. We 
find that the confining potential has a significant effect 
on the DOS in both cases. In the case of quartic confin- 
ing potential (Fig. 1(a)), the divergence in DOS at the 
band edges is found to be suppressed and eventually de- 
stroyed with increasing strength of the potential. Further, 
the quartic potential spreads the DOS over a wide en- 
ergy scale. It is also to be noted that for small potential 
strengths shown in the figure, the low energy part of the 
DOS continues to show similarity to the case of the DOS 
of lattice bosons without confinement. In comparison, in 
the case of harmonic potential (Fig. 1(b)), divergence of 
the DOS at the lower edge of the band is destroyed while 
the one at higher edge is significantly suppressed even for 
a weak harmonic potential, and with increasing strength 
the DOS is flattened to a wide energy scale[23J. In Fig. 
2, we have shown our results on the DOS of a boson in 



confining potentials spread the DOS compared to the pure 
lattice case. Fig. 3 shows the DOS of a boson in a cubic lat- 
tice in a quartic potential. Increasing the quartic potential 
strength clearly has a strong effect on the DOS. Similar 
results are obtained for the harmonic potential case as 
well, as shown in Fig. 4. The dotted lines in Fig. 4 are the 
single particle DOS for an infinite dimensional hypercubic 



lattice whose DOS is [21]: p(E) = exp[-E 2 /(2t 2 )]/^2nt 2 . 
On comparing the changes brought about by the harmonic 
potential in Id, 2d, and 3d, we notice an approximate di- 
mensionality crossover in the DOS for small k. For a Id 
lattice with harmonic potential, the DOS has a finite value 
at the lower band edge and a weak singularity well inside 
the band, which are characteristics of the DOS of a 2d lat- 
tice in the absence of confining potential. For a 2d lattice 
with harmonic potential, the DOS almost vanishes at the 
lower band edge and has a flat region in the middle part 
of the band, which are characteristics of the DOS of a 3d 
lattice in the absence of any confining potential. Finally, 
the DOS of a 3d lattice with harmonic potential is found 
to be close to that of an infinite dimensional hypercubic 
lattice. Hence in DOS, the dimensionality crossovers seen 
are: Id {k ^ 0) -> 2d (k = 0), 2d (k ^ 0) -> 3d (fc = 0), and 
3d (fc 0) — * ood (k = 0). In comparison, such dimension- 



4 R. Ramakumar and A. N. Das: Lattice bosons in quartic confinement 

ality crossovers are not found in the case of lattice bosons 1(a). At low temperatures the bosons are in low energy 

in quartic traps. states. The DOS plots implies that change in the shape of 

the band bottom is not significant for small q. Hence for 

3.2 Ground state occupancy and specific heat small q, the C v curves for different q values are very close 

to each other in the low temperature range. 

In this section we present our results on condensate frac- 
tion, condensation temperature, and specific heat of lat- In Figs. 8-10, we have presented our results for bosons 
tice bosons in quartic potentials. A comparison of the tern- in 2d square lattices in a 2d quartic potential. Compared 
perature dependence of the fractional ground state occu- to the Id case, the condensate fraction is found to in- 
pancy for Id lattice bosons in a Id harmonic trap and in a crease faster with decreasing temperature in the tempera- 
Id quadratic trap is shown in Fig. 5 where we have plotted ture range below Tq. The temperature variation of No/N 
the variation of Nq/N with T/Tq. Here Nq is the boson is seen to be close to the case of pure lattice bosons corn- 
population in the lowest energy level and T is determined pared to lattice bosons in a harmonic trap. This feature 
by setting Nq — and fi = Eq in the number equation is also found in the 3d case discussed later. The results of 
(i.e., by solving TV = J2l=i l/[exp(Ei — Eo)/1cbTo) — 1]). the quartic confinement case is closer to the pure lattice 
We note that the dependence of fractional ground state case since, in the central region of the lattice, the quartic 
occupancy on scaled temperature (T/Tq) is nearly inde- potential is shallower compared to the harmonic confining 
pendent of the strengths of the potentials when lattices are potential. In the low energy range, the DOS of a lattice 
large enough that finite size effects are absent. The lattice with quartic potential is closer to that of a pure lattice 
sizes and the strengths of potentials we have used satisfy compared to the DOS of a lattice with harmonic confin- 
this condition. The dependence of Tq on the strength of the ing potential. Increase in dimensionality leads to smaller 
quartic potential presented in Fig. 6 shows a fast increase number of bosons in the ground state for T > Tq. Fig. 
for small values of q and a monotonic increase for larger 9 shows the dependence of condensation temperature on 
q. The boson number dependence of Tq is found to be lin- the strength of the quartic potential. While the k depen- 
ear (not shown) similar to lattice bosons in a harmonic dence is similar to that found in Id, the magnitude of 
potential [20j. In Fig. 7, we have shown the dependence of T is much lower in 2d compared to that in Id for the 
the specific heat on scaled temperature. Increasing the the same value of q. Fig. 10 shows the dependence of specific 
strength of quartic potential suppress C v except in the low heat on temperature for various strengths of the quartic 
temperature regime. In the low temperature range, the C v trap potential. Unlike in Id, the C v is found to have a 
is nearly independent of the potential strength. This can peak near the condensation temperature. In the interme- 
be qualitatively understood if we regard the DOS in Fig. diate to high temperature range, the specific heat is found 
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to be suppressed with increasing potential strength. Low the overall confining potential, it would be better to use 
temperature C v shows a slight enhancement with increas- a quartic confinement rather than a harmonic one. Fig. 
ing q. When compared with the low temperature C v of 13 contains our results on quartic potential strength de- 
lattice bosons in harmonic traps which has a (T/Tq) 2 de- pendence on the condensation temperature. These results 
pendence, we find that this low temperature part for the are qualitatively similar to Id and 2d, and quantitatively 
quartic trap has approximately a (T/To) 1 ' 7 dependence as similar to 2d than Id. The temperature dependence of C v 
shown in Fig. 11. This can be understood when we con- is shown in Fig. 14 for various strengths of the quartic 
sider the DOS of a lattice boson in quartic and harmonic potential. Increased dimensionality makes the peak at the 
traps. We first recall that the temperature dependence of condensation temperature sharper. Similar to Id and 2d, 
the specific heat of bosons goes as (T/T ) Q for a E^'^ increasing the strength of the potential leads to suppres- 
dependence of the DOS [25 . Now, we find that the low en- sion of C v except in the low temperature region. The low 
ergy part of the DOS of lattice bosons in a quartic trap temperature specific heat has a (T/To) 265 dependence as 
can be fitted with a value of a w 0.8 — 0.9, as a conse- shown in Fig. 11. For lattice bosons in a harmonic trap, the 
quence of which the specific heat exponent is less than low temperature C v has a (T/Tq) 3 dependence. This dif- 
2. For 2d bosons in a 2d harmonic trap, low energy part ference originates from the difference in the DOS for these 
of DOS has a linear E dependence leading to a (T/To) 2 two cases which, in the low temperature range, shows a 
dependence of C v . _£;(<*- !) dependence with a = 3 for lattice bosons in a 

harmonic trap and a ~ 2.5 — 2.6 for the quartic trap. 

Now we discuss our results on 3d lattice bosons in a 
3d quartic trap. The temperature dependence of the con- 

4 Conclusions 

densate fraction, shown in Fig. 12, shows that increased 

dimensionality makes the condensation sharper. The Bose In this paper, we presented results of our calculations of 

condensation is favored by small values of the DOS at the single boson density of energy states, condensate fraction, 

bottom of the energy band. With increasing dimension- condensation temperature, and specific heat of bosons in 

ality, the DOS decreases at and near the bottom of the some one, two, and three dimensional periodic lattices in a 

band and this leads to the sharpness of the condensation, quartic potential. Wherever possible, the results obtained 

Similar to Id and 2d results, the growth of the conden- are compared with the corresponding results for lattice 

sate fraction with temperature is found to lie intermedi- bosons in harmonic traps. In one dimension, we find that 

ate between the lattice bosons in a harmonic trap and the DOS of a lattice boson in a quartic potential continues 

pure lattice bosons. Clearly, if one is looking for phases of to retain its pure lattice form for small potential strengths 

bosons in optical lattices with minimum intervention from unlike the case of harmonic trap. In two dimensions, the 
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Van Hove singularity is suppressed but not eliminated for 5. B. Paredes, A. Widera, V. Murg, O. Mandel, S. Foiling, I. 

small quartic potential strengths in contrast to the effect Cirac, G. V. Shlyapnikov, T. W. Hansen, I. Bfoch, Nature 

of a harmonic potential. In three dimensions, the quartic 429, 277 (2004). 

and harmonic potential are found to destroy flat regions of 6. M. Kohl, H. Moritz, T. Stoferle, C. Schori, T. Essiinger, J. 

the DOS of a pure lattice. In addition, for lattice bosons in Low Tem P- Ph y s - 138 > 635 (2005). 

a harmonic potential, one finds a dimensionality crossover G. K. Campbeil, J. Mun, M. Boyd, P. Medley, A. E. Lean- 
in the DOS, which is not found in the case of a quartic po- hardt > L - G - Maracassa, D. E. Pritchard, W. Ketterle, Science 
tential. The temperature variation of condensate fraction, 313, 649 (2006). 

condensation temperature, and specific heat is found to be 8 ' S - F611in 6' A ' Widera ' T - MMer ' F - Gerbier, I. Bloch, Phys. 

intermediate between those of lattice bosons in a harmonic Rev. Lett. 97, 060403 (2006). 

potential and pure lattice bosons. 9 ' D ' Jaksch < C ' Bmder ' J ' L Cirac ' Q W ' Gardiner, P. Zoller, 

Phys. Rev. Lett. 81, 3108 (1998). 
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Fig. 1. Density Of States (DOS) of a boson in a one- 
dimensional periodic lattice of size 1000 in quartic (a) and 
harmonic (b) potentials. In bottom panel (a): q = 0.25 x 10~ 6 
(dots), 0.25 x 10~ 5 (squares), and 0.25 x 10~ 4 (crosses). In top 
panel (b): k = (dots), 0.0001 (solid), 0.0002 (dashes), and 
0.0003 (dash-dot). All energies are measured in units of t. 
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Fig. 2. The DOS of a boson in a two-dimensional square Fig. 3. The effect of quartic potential strength on the DOS 
lattice of size 150 x 150 in a quartic (a) and harmonic (b) of a boson in a three-dimensional cubic lattice of size 50 x 50 
potentials. In the bottom panel (a): q — 0.25 x 10~ 5 (solid), x 50. 
0.25 x 10~ 4 (dashes), and 0.10 x 10~ 3 (dash-dot). In top panel 
(b): k = 0.01 (solid line), 0.02 (dashed line), and 0.03 (dash- 
dot). 
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Fig. 4. The effect of the strength of harmonic potential on 
DOS of a boson in a three-dimensional cubic lattice of size 70 
x 70 x 70 for the values of k shown on the curves. The dotted 
lines are DOS's for an infinite dimensional lattice in the absence 
of any confining potential, as discussed in the text. 




Fig. 5. Temperature dependence of condensate fraction for 
bosons in Id optical lattices with confining potentials: lattice 
bosons in a quartic trap of strength q — 0.25 x 10~ 6 (dots) 
and lattice bosons in a harmonic trap of strength k = 0.001 
(solid). The lattice size used is 1000 and the number of bosons 
is 600. The values of UbTq for these cases are 3.23, and 6.773, 
respectively. 
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Fig. 6. Quartic potential strength dependence of condensation 
temperature for 600 bosons in a Id periodic lattice of size 1000. 
The lowest value of q in this figure is 0.25 x 10~ 6 . 
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Fig. 7. Temperature dependence of the specific heat of 600 
bosons in a Id lattice of size 1000 in a quartic potential of 
stregnths q = 0.25 x 10~ 6 (solid line), 0.25 x 10~ 5 (dots), and 
0.25 x 10~ 4 (dash-dot). The values of fcsTb for these cases are 
3.23, 6.91, and 14.98, respectively. 
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T/T 

Fig. 8. Comparison of the growth of condensate fraction of 
lattice bosons in 2d: bosons in a 2d square lattice of size 100 
x 100 (dash-dot), lattice bosons in a quartic trap of strength 
q = 0.25 x 10~ 5 (dotted), and lattice bosons in a harmonic 
trap of strength k = 0.01 (solid). The lattice sizes used are 
150 x 150 (harmonic trap case) and 100 x 100 (quartic trap 
case), and number of bosons is 600. The values of ksTo for 
these cases are 0.18, 1.187, and 3.19, respectively. 




q 

Fig. 9. Quartic potential strength dependence of condensation 
temperature for 600 bosons in a 2d square lattice of size 100 
x 100. 
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Fig. 10. Temeperature dependence of specific heat of 600 
bosons in a 2d lattice of size 100 x 100 in a quartic potential 
of stregnths q = 5 x 10~ 5 (solid line), 0.10 x 10~ 4 (dotted), 
and 0.25 x 10~ 4 (dash-dot). The values of fcsTb for these cases 
are 1.187, 1.817, and 2.421, respectivley. 
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Fig. 11. Fits to low temperature specific heats. 2d (bottom 
panel): the dotted line is a plot of 2.6 x (T/To) 1 ' 7 . 3d (top 
panel): the dotted line is a plot of 7.5 x (T/T ( ,) 2 ' 65 
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Fig. 12. Comparison of the growth of condensate fraction Fig. 14. Temeperature dependence of specific heat of 2500 

of lattice bosons in 3d: pure lattice bosons (dash-dot), lattice bosons in a 3d lattice of size 70 x 70 x 70 in a quartic potential 

bosons in a quartic trap of strength q — 10~ 5 (dotted), and of stregnths: q — 10~ 5 (solid line), 10~ 4 (dots), and 10~ 3 

and lattice bosons in a harmonic trap of strength k — 0.01 (dash-dot). The values of fcsTo for these cases are 1.23, 2.48, 

(solid). The lattice size used is 70 x 70 x 70 and number of and 5.63, respectivley. 
bosons is 2500. The values of fesTo for these cases are 0.285, 
1.23, and 2.08, respectively. 
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Fig. 13. Quartic potential strength dependence of condensa- 
tion temperature for 2500 bosons in a 3d square lattice of size 
70 x 70 x 70. 



